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1 Introduction 

There are two major problems in the QCD analysis of polarized inclusive deep 
inelastic scattering (DIS): i) the absence of neutrino data makes it impossible, in 
principle, to determine the non-strange polarized sea-quark densities Au{x, Q^) and 
Ad{x,Q^), ii) the separate determination of the polarized strange quark density 
As{x, Q^) and the polarized gluon density AG{x, Q"^) relies heavily on the QCD 
evolution in Q^ and use of the flavour S'[/(3)i7'-invariance relation 

/ dx [Am + Au + Ad + Ad- 2 As - 2 As] =3F -D. (1) 

Jo 

The absence of a long lever arm in Q^, in the polarized case, and doubt concerning 
the reliability of S'?7(3)i;'-invariance for hyperon /3-decay means that As and AG 
are still rather poorly known [jl|, despite the dramatic improvement in the quality 
of the data in the past few years. 

The direct resolution of these problems must await a series of new machine 
development projects, based on very high intensity neutrino beams, which are most 
unlikely to come into operation before the year 2015. 

In the meantime there is currently a major experimental effort at CERN 0], 
HERA [^ and Jefferson Lab. to study semi-inclusive polarized DIS reactions 

-^ + 1^ ^e + h + X, (2) 

involving the detection of the produced hadron h. 

The theoretical structure for the analysis of such reactions, in both leading order 
QCD (LO) and next-to-leading order (NLO), exists 0]. However we are critical of 
the type of LO analysis carried out thus far by the experimental groups. 

The analysis of semi-inclusive DIS [^ 0] involves both parton densities and 
fragmentation functions. In the LO treatments referred to above, the fragmentation 
functions (FF's) are treated as known quantities from inclusive e+e~ —>■ hadrons, 
and are used in constructing an auxiliary quantity called "purity". However, it 
is well known that in e^e~ -^ hadrons, both in LO and NLO, first, only the 
combinations D^ + D^ can be determined while for the analyses of semi-inclusive 
DIS both D^ and D^ are needed separately, and second, the existing analysis are 
rather ambiguous: a detailed study in p makes a 31 parameter fit to the data, and 
no errors are quoted, and, in a more recent study [0], the FF's differ significantly 



from those in |^, by 40% or more in some regions of z. Under these circumstances 
it is unreasonable to pretend to have an absolute knowledge of the fragmentation 
functions. 

Two NLO analyses based upon a global analysis of the inclusive and semi- 
inclusive data have been attempted 0, f^. In the more recent analysis p the 
authors relax the equality of A'u(x) and A(i(x) and find a preference for a positive 
A'it(x), but effectively no constraint on the sign of !S.d{x). However, in both these 
analyses it is again assumed that the FF's are known exactly: those of |]^ being 
used in [j^, and those of @] in p. 

In addition, in the above mentioned analysis, some simplifying assumptions are 
made about relations between various polarized parton densities. In the following, 
except where expressly indicated, we make no assumptions at all concerning the 
polarized or unpolarized parton densities. Indeed there are persuasive arguments 
from the large-iVc limit of QCD that a significant difference should exist between 
A'u(x) and I\d[x) with \IS.u — A(i| > \u — d\^ and it has been argued that such 
a situation is compatible with all present day data [|l^. Further, bearing in mind 
recent arguments |TT]] that s[x) ^ s{x) and As(x) 7^ As(x), we even refrain from 



the very common assumption of the equality of these densities. 

However, from experience gained in the analysis of inclusive polarized DIS, it 
appears that the parameter space is sufficiently complicated to be able to produce 
biases in the y^ analysis, which can lead to unphysical results. We thus believe 
it to be dangerous in either LO or NLO QCD, to put together all inclusive and 
semi-inclusive data in one global analysis. Rather, what is required, is a working 
strategy, making optimal use of selected parts of the data. 

The aim of this paper is precisely to provide a strategy, in both LO and NLO 
QCD, for the analysis of the semi-inclusive data. We proceed as follows: Informa- 
tion about the FF's is obtained from unpolarized semi-inclusive data. Then this 
information is used to determine the polarized parton densities from polarized semi- 
inclusive DIS. Appropriate use of the information from DIS and e"'"e~ —^ hadrons 
is made as well. We suggest, for example, that one should use as input not just 
a knowledge of the unpolarized parton densities and their errors, but rather the 
polarized isotopic combination 

Aq3{x,Q^) = {Au + Au)-{Ad + Ad) (3) 

which is by now very well determined from polarized inclusive data, and which is 



free from any influence of As and AG. Inclusive e^e" — > hadrons and DIS are used 
also to considerably simplify the NLO expressions for the cross sections of @). 

Further, in stead of dealing with each parton density and FF separately, as is 
often done we work with their singlet and non-singlet combinations. As both the 
parton densities and the FF's enter the semi-inclusive cross sections, we single out 
such observables that are singlet (non-singlet) combinations on both the parton 
densities and the FF's. 

This leads to the fact that we often consider linear combinations of experimen- 
tally measured quantities, and it may be objected that thereby we are dealing with 
experimental observables with possibly large errors. It has to be understood that 
that is a reflection of the true situation and not an artifact of our approach. To take 
an absolutely trivial example, suppose Ei and E2 are two experimentally measured 
functions used in an attempt to determine the theoretical functions Ti and T2, where 
El = Ti + T2 and E2 = Ti ~ T2. Now, if it happens that Ei ~ E2 and if we write 
^1,2 = (1/2) (-El ±-^2) then T2 will be very poorly determined. This is unavoidable. 
It does not help to do a best fit to Ei and E2 with some parametrisations of Ti and 
T2. Inherently the result for T2 will have a large relative error. 

Thus we believe that any relatively large errors occurring in our manipulation 
of the experimental quantities reflects a genuine and unavoidable imprecision in the 
determination of certain theoretical quantities. 

In order to minimize systematic errors experimentalists prefer to consider asym- 
metries or ratios of cross sections where the detection efficiencies should roughly 
cancel out, e.g. ratios of polarized to unpolarized DIS or ratios of polarized to 
unpolarized semi-inclusive for a given detected hadron. We appreciate that this 
is a fact of experimental life. However we wish to stress that a large amount of 
information is lost in restricting oneself to only these ratios. It is vitally important 
to gain control of the systematic errors in detection efficiencies, and although we 
try as far as possible to deal with the favoured kind of ratios we will be forced also 
to consider other types of cross-section ratios. 

Throughout the rest of this paper we assume that a kinematic separation is 
possible between hadrons produced in the current fragmentation region and those 
produced from the target remnants. We consider only the current fragmentation 
region so that our formulae apply only to this region and fracture functions [|1^ 
play no role in our discussion. 

It has to be stressed that there is a huge difference in complexity between the 



LO and NLO treatments. Thus it makes sense to utilize the LO approach, provided 
appropriate checks, (which we suggest) are carried out. 

Our analysis proceeds in a step-by-step fashion. Firstly we describe a generic 
test for the reliability of a LO treatment. Assuming this to be successful we present 
the analysis in LO. In LO the information on the polarized valence quark densities 
Awy and Ady and on the breaking of SU(2) invariance of the polarized sea is 
obtained without any knowledge of the FF's. However, it should be clear that any 
information on the sea quark densities in LO cannot be reliable, since they are 
expected to be small and thus comparable to the NLO corrections. We consider 
also the experimentally difficult case of (f) production, which seems to be the best 
way to get an accurate determination of As + As in LO. We discuss also how, in 
principle, one can test s{x) = s{x) and/or As(a;) = As(x). 

In the NLO treatment we show in a new way how information on e~^e^ —>■ h + X 
and inclusive DIS can be incorporated directly so as to considerably simplify the 
NLO expressions for semi-inclusive cross sections. This key result is presented in eqs. 
( |69D and (|70D . Next we discuss the fragmentation combination D'^ — D'^ and use 
it to evaluate Awy and /S.dv in NLO. Then we obtain expressions for _DJ +-0^ , for 
D"^ and for D^^ +D1'' . With these we are able to obtain (Am— A J) in NLO. Finally 
we have a set of 2 equations involving 3 unknown functions, (Au + An + Arf + Ad), 
(As + As) and AG. An accurate determination of all three functions would require 
data over a presently impossibly wide range of Q^. We thus suggest using here for 
AG its determination from charm production [|l^]. It should then finally be possible 
to get an accurate assessment of As + As in NLO. Lastly we consider the evaluation 
of s(x) — s(x) and As(x) — As(x). 

2 A strategy for semi-inclusive DIS 

In NLO QCD, the expressions for semi-inclusive cross sections involve convolutions 
of parton densities and fragmentation functions with (known) Wilson coefficients . 
Our lack of knowledge of the errors on the FF's will make it difficult to assess the 
accuracy of the parton densities which we are trying to determine. In the LO QCD 
approximation, on the other hand there are no convolutions, but simple products 
only (independent fragmentation), and it becomes possible to construct measurable 



combinations of cross sections in which the FF's completely cancel out |jT5|, piq] . 



However it is not clear how reliable the LO is. 

We believe it is quite safe when determining the large Am and Ac? densities, but 
could be quite misleading for Au, Ad and As. In any event it is absolutely essential 
to test independent fragmentation in order to have a feeling for the errors on parton 
densities obtained via the LO formalism. 

In LO, the structure of the expressions is generally of the form 

parton density Aq{x, Q^) = experimental observable E{x, z, Q^) (4) 

or 

fragmentation function D(z,Q'^) = experimental observable E(x, z,Q'^). (5) 

In both cases the characteristic feature of the LO treatment is that the RH sides, 
which can in principle depend on (x, z, Q^), should only depend on two of these, 
either {x, Q'^) or (2;, Q"^) respectively, so that there is an independence of the third 
variable, which we shall call the passive variable. 

Every expression of the form (^ or (|5|) should be tested for dependence on 
the passive variable. If a significant dependence is found it does not mean that 
the LO analysis must be abandoned, but it suggests that the variation with the 
passive variable be used as an estimate of the theoretical errors, 5TH[Aq{x,Q'^)], 
(5th[-D(2;, Q^)] on the sought for quantities. 

In Section 7 we discuss a strategy for the analysis in NLO. 

3 Parton densities in LO QCD 

It is useful to introduce the following notation for semi- inclusive processes: 



and 



Atio^ \l + {l-yyj dxdydz 



^,.^<p + 'f 



Aira"^ \2 — y 



d^a^^ d^ 



o'l 



dx dy dz dx dy dz 



(7) 



where P^ and l^ are the nucleon and lepton four momenta, and cta/^ refer to a lepton 
of helicity A and a nucleon of helicity ^. The variables x, y, z are the usual DIS 
kinematic variables. 



Then in LO the cross sections for the semi-inclusive production of a hadron h 
have the simple y-independent form 






(9) 



<i,<i 



where the sum is over quarks and aniquarks, and where D'l is the fragmentation 
function for quark or antiquark i to produce h. 

We consider sum and difference cross sections for producing h and its charge 
conjugate h on both protons and neutrons, and define 



AAi%\x,z,Q') 
AA';tl{x,z,Q') 



Aa^ ± Aa^ Aa^^^ 

p,n ^ p,n '—^^p,n 



p,n P,n 



p,n 



;^h±h 



p ^ n 

pr-h±h _|_ pj-h±h 
^p ^n 



(10) 



111 



For inclusive unpolarized and polarized DIS cross sections we use the notation: 

2y' 



^DIS^^P + l) 



2 / o„.2 \ ^2^DIS 



and 



Aa 



~DIS 



Ana^ \1 + {1 — yy J dxdy 
x{P + lf f y 



d^a^'J d^a^'^ 



dx dy dx dy 



47ra^ \2 — y ^ 
Then in LO we have the ^/-independent expressions: 

Aa^^^(x,Q2) = 2^f(x,Q2)Uo = 5^ejAg,(x,Q2 



(12) 
(13) 

(14) 
(15) 



9,9 



In addition to ([ToD and ([Tl|) we consider the ratios of sum and difference hadron 
yields for the unpolarized semi-inclusive and inclusive processes: 



%-h _i_ x-h 



^p,n l^! ^5 W ) ~ ~nT.<^ ' -^pztnl-^'^'W J 






~h±h _j_ ~h±h 
p ^ n 

Z-DIS -U ^DIS ■ 
^p ^ n 



(16) 



(It is equally good to use a sum over any set of hadrons h and their charge conjugate 
h.) 



3.1 Testing LO QCD 

Using only charge conjugation invariance it is easy to show that [|16 

This is a key relation for testing the reliability of the LO. The RHS is completely 
known from inclusive DIS and is independent of z. The LHS, in principle, depends 
also upon z and upon the hadron h. Only in LO (or in the simple parton model) 
should it be independent of z and of h. It is thus crucial to test this feature. 

To help with statistics it is also possible to formulate an integrated version of 
(|17|). This is given in |]T6|. For the rest of this section we assume that the test ( |T^ ) 



has been successful and proceed with the analysis in LO. 

3.2 The valence quark densities in LO 

The polarized valence quark densities can be obtained from tt^ production, assum- 
ing only isospin invariance: 






15 
Mv = 



15 

For the case of K^ or A, A production, if one makes also the conventional as- 
sumption that s = s and As = As one has in addition: 

Ady = l{K + t/v)AAf+V^"-K-rfv)A<_V^'}. (19) 

Note that we have not assumed D^ = Df , although that is suggested by the 
absence oi a d quark in the leading Fock state of K"^. Indeed the above equality 
can be tested (see Section 5.2). 

For isoscalar hadrons like A, A again assuming s = s and As = As, one finds 

Auv = ^{A{Auv + dv)AA^-~^-{Adv + uv)AAy^^ 

^ |4(4dy + uv)AAy'' - {Auv + dv)AAy^^ . (20) 



Adv ■ 

^ 15 I 
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We shall comment in Section 5.5 on the situation if one does not assume s = s 
and As = As, where we suggest a method for estimating if the failure of these 
equalities is serious or not. In any event, the safe way to obtain Auy and Ady is 
via TT^ production, eq. (|T^). 

Once again the reliability of these LO equations can be tested by checking that 
the RH sides of (p!8|), (|19|) and ( pOf ) do not depend on z. If it is found that for a 
given x-bin the RH sides vary with z by some amount 6th[^uv], STni^dv]: then 
these could be regarded as an estimate of the theoretical error at this x value. 

4 Use of Aqs{x, Q'^) from polarized inclusive DIS 

At NLO the spin dependent structure functions for protons and neutrons are given 
by: 



q=u,d,s 



{Aq + Aq)®[l + ^^6C,] + 



2^ 



+ ^^AG®6Cg 

In 



(21) 



involving a convolution of polarized parton densities with known Wilson coefficients. 
For the neutron, (?" is obtained by the replacement 

An ^^ Ad (22) 

p 

Now it is clear that one can obtain information only on the combinations: 

Au + Au, Ad + Ad, As + As, AG. (23) 

and that it is impossible to obtain separate information on the valence and non- 
strange sea quark polarizations from inclusive, neutral current, polarized DIS. 

In our semi-inclusive analysis we shall use as a known qauntity only the non- 
singlet combination of the polarized parton densities Aq^, eq.(|^). Now that there 

^However, sometimes it is convenient to parametrize separately the valence and sea quarks and 
to assume, for example, Au = Ad = AAs, where A is a free parameter. It should be obvious then, 
that any claim that the x^ analysis favors some particular value of A must be fictitious and a 
consequence of some hidden bias in the minimization procedure. Yet such claims have been made. 



is such an improvement in the quahty of the neutron data we beheve this quantity 
is very well constrained directly by the inclusive data. For one has, from (0) 

gl{x, Q') - <(x, Q') = ^ Ags ® (l + ^^^C,^ (24) 

and Aq^{x, Q^) is determined without any influence from the less well known quan- 
tities As and AG, either in LO or in NLO. Of course if the semi-inclusive analysis 
is done in LO one must use Ag3(x, Q^)|lo. Note that we do not use information 
on (Am + Am) or {Ad + Ad) from inclusive DIS, since these are subject to the less 
known strange quark and gluon effects. 



4.1 SU(2) symmetry of the sea quark densities in LO 

One has 

[Auix, Q') - AJ(x, Q^)]lo = \ [Ag3(x, Q^) + Arfy(x, Q^) - Auv{x, Q^)]^^ . (25) 

Eq. (p5D determines the SU(2) symmetry breaking of the polarized sea without 
requiring any knowledge of the unknown Ag and AG. A possible test for SU(2) 
breaking for the polarized sea densities that does not require any knowledge of the 
polarized densities was given in |jT6l . 



In order to determine the polarized sea quark densities separately we need one 
more relation, namely the value of 



For then 



Aq+{x,Q^) = Au + Au + Ad + Ad. (26) 



(Am + Ad)Lo = ^{Aq+ - Auv - Adv)Lo (27) 



which combined with ( [25| ) yields Au and Ad separately. 

Although each term on the RHS of (pSf ) and ( pT]) should be well determined 
in LO, the corresponding linear combinations are expected to be small and may 
thus be very sensitive to NLO corrections. An indication of the sensitivity may be 
inferred from the fact that in inclusive polarized DIS, As{x, Q"^) changes by roughly 
a factor of 2 in going from LO to NLO or when one changes factorisation schemes 
from Ms to AB or JET. 



Note that determining, say, Am via Am = |[(Am + An) — Auy] is unreliable 
since determination of (Aw + Au) from inclusive DIS requires a knowledge of As. 

In order to determine Ag+(a;, Q^) it will first be necessary to extract some in- 
formation on the FF's. 

4.2 Fragmentation functions in LO 

1. From measurements of the ratios Rpln oi the semi- inclusive to inclusive DIS 
cross sections on protons and neutrons for any given hadron h, it is feasible in 
LO to learn a great deal about the FF's D'^ + D^ = D^^^. This combination is 
measured also in e~^e~ -^ hadrons. 

Analogous to the polarized case, we define 

q-i{x,Q^) = u{x,Q^) + u{x,Q^)~d{x,Q^)-d{x,Q^) (28) 

q+{x,Q^) = u{x,Q^) +u{x,Q^) +d{x,Q^) +d{x,Q^). (29) 

which are well determined from inclusive DIS data and which thus can be taken as 
known quantities in the semi-inclusive analysis. 



Using data on unpolarized semi-inclusive DIS we have, in LO, 

'4D'^-'\z,Q')-D'/\z,Q')]. (30) 



9 
ph+h ^ f_ 

-'^p—n q 



• For the case of pions, kaons and A, when SU(2) invariance can be used this 
simplifies to 

i?;:r" = 2Df+^-{z,Q') = 2Df+--{z,Q'). (31) 

• For K mesons and A hyperons, 

Rp^n''' = 2D^'^^\z, Q^) = 2Df^^\z^ Q2), (32) 

where the superscript K stands for the sum over all produced kaons: 

K = K+ + K- + K^ + K^. (33) 

It will also be of great interest to compare these FF's with those obtained from 
e^e~ — > hadrons [^ ^ and those used in Monte Carlo models. 

10 



2. Given that My(x, Q^) and dyiXyQ"^) are well determined from inclusive DIS 
one can proceed further to obtain the other combinations of FF's D^ — D^ = D^^^. 

• One finds for vr^ 



r^'T^ —TV 



9< 



4:Uv — dy 



Auv — dv 



18 (Ff 



1 DTT'"— 7r 



Adv — Uy 4:dv — Uy 

Combined with (|31|) we have expressions for D„ and D^ separately. 
For K^ one obtains 



(34) 



4D 



K+-K- 



D 



K+-K- 



Q{a 



~K+-K- 



^ n 



Uy — dy 



1^ 



JDK+-K- 
'LO^n 



, (35) 

Uy — dy 

It is usually assumed, and this is presumably a very good approximation, that 



D 



K+ 



D 



K- 



(36) 



in which case ( pSf ) can be read as an expression for D^ ^^ . It is not possible 
to test relation (|3^) without taking s = s and/or As = As, but such an 
approach is hard to justify given that any failure of (PUf ) is presumably very 
small. 



For isoscalar hadrons like A, A 



D 



A-A 



D 



A-A 



(-^1 Jlo ^P 



A+-A- 



v-^1 )lO ^n 



A+-A- 



Uy — d 



V 



(37) 



Given that Auy and Ady are determined via (|TB]) we can write analogous expression 
for the above D^ functions using the polarized data. 
Of course all the expressions (|30D 



S), (H), (H, (H) and (13), being LO 
results, must be tested by demonstrating that the RH sides are essentially indepen- 
dent of the passive variable x. 

Now that we have determined DZ ^'^ we can determine DZ ^'^ in LO via 



DZ 



^{Ff + F^)^oR. 



L O ^ ^p+n 

2 (s + s) 



- 5g+/^r 



(3^ 
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Similarly, since D^'^^ is determined via (|32D, we can find D^'^^ from 

4.3 The non-strange sea quark densities revisited, in LO 

Now that we have determined Z^^ ^'^ and DJ '^'^ in LO we can, in principle, de- 
termine Ag+(x, Q^) and As(x, Q'^) + As{x, Q"^) from the semi-inclusive and inclusive 
relations, in LO, 

.f + .r ^ ^^-^--^^3^^^^^'^ (40) 

.,.++.- ^ 5Ag+Df+--+2(As + As)Z}f+-- 

P+n 5g+D-++-- + 2 (s + s)Df +-- ■ ^ ' 

Such a determination of Ag+ in LO is likely to be reliable, but As + As and 
the individual A-u and Ac? obtained in LO via (^) may be subject to significant 
uncertainty. 

We note that there is an alternative way to determine Ag+, but it requires the 
ability to detect K^ . In that case one has, in LO, 

. _ AOp + A(Tn ^_^ , . 

^?+ — 1+ ^K++K-~(KO+KO) , ~K++K--(KO+KO) ' v^^'' 

0-p + (Tn 

4.4 The strange quark density As + As in LO 

The approach to As + As in Section 5.3 and the approach discussed in |I6[ are 
unlikely to be reliable. The problem is that for production of pions the strange quark 
contribution is "doubly small", since e.g. one must compare AmD^ with AsDJ in 
which both |As + As| <^ |An| and \D1\ <^ \D'^\. For kaons and A hyperons it is 
somewhat better in that \D^\ k. \Df\ and \D^^~^\ = \D^^^\ ^ \DJ+^\, but this is 
similar to the situation in inclusive DIS where we know that the LO determination 
of As + As is quite unreliable. 

The only possibility we can see for a reasonable determination of As in LO is 
via production. For in this case one has, |As + As| ^ |Aii| but presumably 
\Df\ ^ |D^| so that the strange and non-strange quarks are on equal footing. 
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One has, by charge conjugation invariance Df = Df, and it should be quite safe 



to take Dt = Dl = D^ = D'l^. One then obtains in LO 

As + As _ 3A<Tp+„ - 5 y-^j Acrp_„ 

Moreover one has expressions for the fragmentation functions as weU: 



2(s + s-)H+'^ HI 



Dt = lTTrV-AH^n-^[-]^t-n\ (44) 



Dt = -^^. (45) 

Q3 



As always expressions (^3]) - (^51) must be tested for non- dependence on the relevant 
passive variable. 

Finally we note that the same eqs. (^) - ( ^5D hold for K, A+A and vr-production 
if the superscript (p is replaced by i^, A + A and vr^ + tt^, respectively. And though 
the production rates are higher, the sensitivity to the strange quarks in the i^, A + A 
and TT-productions is lower. 

4.5 Concerning s = s and As = As in LO ^ 

In the analysis of DIS it is conventional to assume s = s and As = As. However 



there are models and arguments ||TT| which suggest that these equalities might not 
hold. 

Within the limitations of the LO we can test these relationships via {K~^,K~) 
and (A, A) production. One has for the unpolarized case, assuming D^ ~^ =0 
several different possibilities: 

{s-s)Dr-''- = i8(Ff),^i?r-^--4«v/^r-^- 



9{uvar-'''-dva^^-^-} 



. (46) 

Uy — dy 

Then for the polarized case, given that Awy and Ady are known from (0) and 

(s — s)Df ~^ is determined, we can proceed to determine (As — As)Df- ~^ 

'^We are grateful to M. Anselmino, M. Boglione and U.D'Alesio for drawing our attention to 
this issue 
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fromAA^^"^ or AA^^~^ : 



A /l-f _ 



. ,i.+-i.- ^ 4AdvDf "-^- + (A. - A.-)Dr-^- 



u 
^K+-K 



D^ ^ is assumed to be known through (|35D. 

For isoscalar hadrons hke A, A we have the possibhties 

{s--s)Dy^ = 18 {Ff),oRy'^- {Any + dv)Dt-'^ 



I {4uv + dv) a^-^ - {My + uy) a^"^} , (4^ 



Uy — d 

where D^~^ is assumed to be determined in (pTl). Here and in eq. p6|) it is 
clear that a^ ~^ and a^~ can be replaced by the corresponding experimentally 
measured quantities 2 {Ff)^^R^^~^' or 2 {F(^)^^R^-^. Then (As - As)D^-^ 
can be determined via AA^'^ or AA^~^: 

^^A-A ^ (4A^.y + Ady)Dt^ + (As - Ars)D^-^ 
^ {Auy + dy)D^-~^ + {s- s)D^~~^ 

.J, {AAdy + Auy)Dt'' + (A3 - As)Dt'' ,,.. 

{4dy + Uy)Dt^ + {s-s)D^-^ ^ ^ 

Although (^), (|i7|), ( ^HD and (|^), being LO expressions, cannot be expected 
to yield accurate values for s — s and As — As, they should nonetheless enable one 
to say whether they are compatible with zero since D^~^ should be relatively large. 
Note that s — s and/or As — As different from zero would break the independence 
of the RH sides of (^) and (pQ]) on the passive variable z. Of particular interest is 
the speculation that As ^ —As but s ^ s, the consistency of which could perhaps 
be tested from (|Up, (g^), (|Hp and (||). 



5 Semi-inclusive analysis in NLO QCD 

The situation in NLO 0, Il7[ is much more complicated than in LO, since factor- 



ization is replaced by convolution, and it is also more complicated than inclusive 
DIS in NLO since here one has to contend with double convolutions of the form 
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q®C®D and Ag ® AC ® D for the unpolarized and polarized cases respectively, 



where C and AC are Wilson coefficients ffist derived in |]l3 and 
The double convolution is defined as 

where the range of integration is given as follows: 



Ii : The range is 



X 



X + {1 — x)z 



< x' < 1 with z < z < 1 



(50) 



(51) 



I2 : In addition to (|5lD there is the range 
X < x' < — 



X 



X 



x)z 



x( 1 — X ] 

with —, '- <z'<l. 

x'[l — x) 



(52) 



Note, that contrary to the case of the usual DIS convolution, the double convolution 
g C -D is not commutative. 

We shall frequently encounter expressions of the form 



a.. 



/vr 



(53) 



corresponding to the LO plus NLO corrections. The fiavour structure of the results 
becomes much more transparent if we adopt the following symbolic notation: 



«« 



qD + —q ®C ®D = q 
/vr 



1 + ®^C® 



D. 



(54) 



Then the semi-inclusive polarized cross section Aa^ defined in (|^) is given by 



AaJ = 5^e2AgJl + ®^AC,,C 



K + 



+ {j2^'^l] ® ^^C,, ® D^ + AG ® ^AC,, ® iTe^D^Jibb) 

where the sum is over quarks and antiquarks of fiavour i and parton densities and 
fragmentation functions are to be taken in NLO. 

For the unpolarized semi-inclusive cross section in NLO it is not possible to 
completely factor out the y-dependence. Consequently a'* defined in @ will depend 
upon y in NLO, in contrast to the LO situation in (0). 
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In the notation of the seminal paper of Graudenz |jT^ the cross section is a 
sum of what he refers to as "metric" (M) and "longitudinal" (L) terms, with corre- 
sponding Wilson coefficients C^ and C^ . A further complication is that the Wilson 
coefficients are different in the two regions of integration Ji and l2- Thus we have 
coefficients C^^\ C^^ with j = 1, 2. 

We then define the y-dependent combinations of Wilson coefficients: 

CJ, = Cl^ + [1 + 47(Z/)] Cl^, (56) 

q, = Ci^' + [1 + 47(1/)] Ci^ (57) 

q. = ^^f + [1 + 47(2/)] Cl^, (58) 
where 

Then the unpolarized semi-inclusive cross section can be written in a form anal- 
ogous to the polarized one: 



~^i = ^^^^^1+®^'^'?'?^ 



K + 

HI 



+ E^'^^ ® ^C,,®D% + G®^C,,® Y.<Dl . (60) 



27r "^ ^ 27r 



In ( |60D we have used the symbolic notation: 

q^ ® ^C,, ® Dj = //^ ® ^^;. ® ^.' + £ ^^ ® S^''^ ® ^'^ ' ^^^^ 

and analogously for C^^ and Cgq. 

Note that in NLO the unpolarized inclusive cross section, a^^^ (|T2|) is given by 

a^^^ = 2Fi [1 + 2-f{y)R] , (62) 

where R is the usual DIS ratio of longitudinal to transverse cross sections. Note that 
strictly speaking here and throughout the rest of this paper R should be replaced 
by R — > i R — ^^3#- ) / ( 1 + ^^^2^ ) , since the correction terms may be important 
for low values of Q^. 

As in the LO discussion we doubt the reliability of a global NLO analysis of 
inclusive and semi-inclusive data, and we suggest that one should feed into the 
semi-inclusive analysis as much reliable information as one can from other sources. 
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As we shall see there is the oft found opposition between what is simple the- 
oretically and what is simple experimentally. However, if the systematic errors in 
detection efficiencies can be brought under control then we can make remarkable 
theoretical simplifications and we can then extract a vast amount of information 
from the semi-inclusive data. This is an important experimental challenge as will 
be seen from the power of the results given below. 

5.1 Simplification of the semi-inclusive NLO results 

Bearing in mind the NLO result (^) for g^ for polarized DIS and the analogous 
result for Fi in unpolarized DIS, we see that in the second term of (|60D we may 
make the replacement, correct to NLO, 

J2eh-^2Ff (63) 

i 

and similarly, in the analogous eq. for AcXp in (|55D : 

J2e^^q^-^2g',. (64) 

i 

Note that here and throughout the rest of this paper Ff '" and gi'^ are the experi- 
mentally measured structure functions. 

Further in the reaction e"'"(p+) + e^(p_) —>■ h + X in the kinematic region when 
we can neglect Z'^-exchange effects, we have 

a''{z, cos^, Q2) = -(1 + cos^ e)a!^{z, Q^) + -Al- cos^ e)al{z, Q^) (65) 

8 4 

where Q^ = {p+ +p_)^. 
In NLO QCD one has 

4(^, Q') = 3ao JE ^.'^i ® (l + g^j) + 



(66) 



where the sum is over quarks and antiquarks, the C^'s are Wilson coefficients, and 






^0 = ^^TW- (6'^) 
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Then, correct to the required NLO accuracy, in the third term in (^0]), and in 
its polarized analogue (|55|) , we may make the replacement 



H-^Di 






(6^ 



where a^ is the experimentally measured cross section. 

Hence, for the unpolarized and polarized semi-inclusive cross sections, in NLO 
accuracy, we have 






E 



ehi 



1 + ®^Cg,c 

/vr 






a. 



+2Ff ® ^C,, 

ZTT 



1 

3ao 



Dg + IT—G Vy TTZ-^gq 



27r 






(69) 



Aa, 



:^h 



5Z^*'^^* 



a. 



1 + ®-^ACg,C 






+2g',' 



2tt' 



AC, 



qg 



D% + 



'G 



3an 



-AG 



27r' 



AC, 



gi 



(Jn 



(70) 



Given that the unpolarized gluon density is reasonably well known, the last term 
in (^9]) can be considered as a known quantity. In the following we take as known 
quantities the NLO values for g+(a;, Q^), 5'3(a;,(5^), My(x,(5^), dy{x,Q'^\ G{x,Q'^) 
and Aq3{x,Q^). 



5.2 The polarized valence densities in NLO 

Using charge conjugation invariance one obtains, for semi-inclusive pion production 



DTT ' — TT 






18Ff[l + 27(y)i?P] 

[Ady - uv][l + ®{aj2n)Cgg®]Df-^' 

18Ff[l + 2^{y)Rp] 



(71) 



The only unknown function in these expressions is D'^ ~'" [z^Q"^)^ which evolves 
as a non-singlet and does not mix with other FF's. A x^ analysis of either or 
both of the equations (|7T1) should thus determine D^ ~'^ in NLO without serious 
ambiguity. 
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Assuming D,^ '" is now known, one can then determine A-uy and Ady in NLO 
via the equations 



^^.+-.- ^ (4A..V - Arfv)[l + 0(«,/27r)AC,,0]Df -- 



r+- 



where, of course, Awy and Adv' evolve as non-singlets and do not mix with other 
densities. Eqs. (^) and ( [T5| ) determine the densities iSxiy and A(iy in NLO without 
any assumptions about the less known polarized gluon and sea densities. 

5.3 SU(2) symmetry of the sea quark densities in NLO 

Once Awy and Arfy are known in NLO we can calculate 

[Au(x, Q2) _ aJ(x, Q'')\nlo = ^[^qsix, Q2) + Adv{x, Q^) - /\uv{x, Q^)]nlo. (74) 



Eq. ([7^) determines the breaking of SU(2) symmetry for the polarized sea densities 
in NLO without requiring any knowledge of Ag and AG. It will be interesting to 
compare the values obtained from (|7^ with information on Au{x) and Ad{x) which 
will emerge from Drell-Yan and W^ production experiments at RHIC [|18|. 

The separate determination of Am and Ad requires knowledge of Ag_|_, defined 
in (0), in NLO. 

Note that determining, say, A-u via A-u = |[(A'u + A-u) — A-uy] is unreliable, 
since the determination of (Am + An) from inclusive DIS involves knowledge of 
gluon and strange quark densities. 

As in the LO case we can only determine Ag+ after obtaining some information 
about the fragmentation functions. 

5.4 Fragmentation functions in NLO 

We consider the sum for the unpolarized production of tt^ and tt^. We have 



6 [Ff (1 + 2-f{y)Rp) -FI'{1 + 27(1/)/?" 



K n ' = ^ ,^r>r TTTTl FT^^^ TTTTT^ ' (75) 
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aa;!+- = *7 r^ ■ (76) 

ga { [1 + ®i^C„®] Dr+-" + ®tC,, ® Dl^-' - 

The only unknown functions in these relations are DJ ^'^ and D^ ^'^ , which will 
mix with each other under evolution. Thus it should be possible to obtain them 
from a x^ analysis of (^) and (|76D. 

Once we know DJ[ ^'^ and utilise D"^ "'^ from Section 6.2 we clearly have 
access to both Z)J[ and D'^ . 

Note that if K^ can be detected one can obtain information on D^ and Dq, 
where K = K^ + K~ + K^ + K^. One simply replaces the labels vr"'' + vr^ by i^ 
everywhere in (l75| ) and (|76D . Analogous equations hold also if 7r+ + vr" is replaced 
by A + A, if A is detected. 

Returning to the case of 7r+ + tt^, the ratio -Rp+^'^ allows the determination of 
the only unknown function D^ ^'^ . We have 

= { (5g+ [1 + ®tC„®] Df+-~ + 2{s + s) [l + ®fjC„®] Df +-') + 

+18(Ff + F^) ® fjC„ ® DS'^^" + ^G ® §C„ ® a?"+-- } / 

/6 [Ff (1 + 2^{y)RP) + Fi" (1 + 2^{y)R^)] (77) 

Under evolution Z}J[ +'^ and Z^J ^'^ mix with D^ ^'^ , but this is not a problem 
since the latter is supposed to be known. It would be interesting to compare these 
results with those from e^e~ — >■ hadrons, obtained recently in ^. 

Again analogous sets of equations holds for kaon and A, A production. One 
simply replaces 7r+ + tt^ by i^ or A + A and this allows the determination of the 
only unknown function D^ or D^^^, respectively. 

5.5 The sea-quark densities in NLO 

With the NLO knowledge of D^ ^'^ , Z^J ^'^ and Dq ^'^ we are now in a position 
to try to determine Ag+(x, Q^), As{x,Q'^) + As{x,Q^) and AG{x,Q'^). We have, 
in NLO, 



9' + 9i=T^ I5A9+ + 2(As + As)] S 1 + ^f^SC, + ^^AG ® SCrj, (78) 
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and 



A ~7r++7r _|_ A ~7r++7r 

A<r" = \^ I (79) 



where 



5Ag+ [1 + ®f^ACgg®] DZ^+^' + 2{As + As) [l + ®f AC,,®] Df+^' 
+18 K + 9?) ® ^AC*,, ® D5^+-" + A AG ® ^AC,, ® a?^+-" (80) 



and 






(5g+ [1 + ®tC,,®] Df +-- + 2(s + s) [1 + ®tC,,®] Df' 



+ 18 (Ff + Fi") ® ^C,, ® DS++-" + ^G ® ^C,, ® aj"'+'^". (81) 

Note that an analogous set of equations hold for 7r+ + vr^ replaced by i^ or A + A. 
Eqs. (|78|) to (|8l| ) contain three unknown functions Ag+, As + As and AG, which 



nonetheless can all be determined in principle because of their different evolution 
in Q^. However, to be at all efficacious such a determination would require a huge 
range of Q^, far larger than is available in present day polarised DIS. 

On the other hand there is a direct and superior method for obtaining AG, 
namely via cc production. This is one of the major goals of the COMPASS exper- 
iment at CERN. We shall thus assume that AG has been determined, so that the 
last term on the RHS of (^) may be taken to be known. 

It should be then straightforward to determine Ag+ and As + As in NLO from 
a x^ analysis of (|78D to (|8lD in which the evolution of Aq^ and As + As would 
involve mixing with the supposed known AG. 

An independent determination of Ag+ and As + As could be obtained by com- 
bining (^) with AAp_^,^ for kaons and Ay4p^„ for A, A production. 



Once Ag_|_ is known in NLO, we can obtain the individual An and Ad from (173) 



5.6 Concerning s = s and As = As in NLO 

It is possible to get some information on s — s and As — As in NLO. 
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We have, assuming D^ = D^ 



K+ _ nK- 



Rp 



K+-K- 



s 



4uv [1 + ®tC,,®] D^^-^- + {s-s)[l + ®^C,,®] Di 

/18Ff [1 + 2j{y)RP] (82) 



= [Ady [I + ®tC,,®] D^^'^- + {s-s)[l + ®^C,,®\ /^r-^"} / 

/18F- [1 + 27(|/)i?"] (83) 

These two equations, taken together with those for AA^ ~^ and AA^ "^ : 

= (4Amv [1 + ®tAC,,®] D^^-^- + (As - As) [l + ®§AC,,®] Df^'^'^ / 
I {^uy [1 + ®^C,,®\ D^^-^- + {s~s)[l + ®^C,,®\ Df^'^-) (84) 

= (4Ady [1 + ®f^AC,,®] D^^-^- + (As - As) [l + ®i^AC,,®] D^^-^'^ / 
I {^Uy [1 + ®tC,,®] dT'""' + (s - s) [1 + ®fC,,®\ Z^r-^-) (85) 

provide four equations for the three unknown functions i^s — s) ® Df ^^ , (As — 
As) (g) Df ~^ and D^ ~^ , so that, in principle, all can be determined via a x^ 
analysis. 

In addition one has 

pA-A _ 
Kp - 

= {{4uy + dy) [I + ®tC,,®] D^^ + (s - s) [1 + ®fSm®\ Dt~^] I 

/18Ff [1 + 2^{y)RP] (86) 

and 

pA-A _ 

= {{Ady + Uy) [I + ®tC,,®] Dt'"" + (S - S) [1 + ®§C,,®] D^^} / 

/18Fi" [1 + 27(|/)i?"] (87) 
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together with AA^'^ and AA^~^: 

AA^'~^ 



p 



= {{4Auv + Adv) [I + ®t AC,,®] Dt~^ + (As - As) [l + ®f^AC,,®] D^^) / 
/ ((4ny + dv) [I + ®i^C,,®] Dt'^ + (s - s) [1 + ®f^C,,®] D^^) (88) 

AA^^ = 
= {{4Ady + Ady) [I + ®^AC,,®\ Dt^ + (As - As) [l + ®^AC,,®\ 0^'"") I 

I ((4rfy + dy) [1 + ®i^C,,®] Dt^ + (s - s) [1 + ®tC,,®] Z),^-^) (89) 

These provide four more equations but only two new unknown functions D^~^ 
and D^~^. The system of eight equations (0) - (^) therefore over- const rains 
the unknown functions and might, hopefully, allow a reasonable determination of 
the relation between s — s and As — As and whether or not s{x) = s{x) and/or 
As(x) = As(x). The actual determination of s — s or As — As requires knowledge 
of either D^~^ or Df ~^ . These could be taken from the study of e^e^ -^ hX. 
This completes the determination of all the polarized densities in NLO. 

6 Conclusions 

We have argued that the present LO QCD method of analysing polarized semi- 
inclusive DIS, using the concept of purity, is quite unjustified. We have also argued 
against attempts at a global analysis, either in LO or in NLO QCD, based on the 
combined data on polarized inclusive and semi-inclusive DIS and taking as known 
exactly the relevant FF's. 

Instead, we have presented a strategy for a step by step evaluation of the po- 
larized parton densities and fragmentation functions from semi-inclusive data using 
selectively chosen information from inclusive DIS reactions. 

In our approach the usually made simplifying assumptions about relations be- 
tween Am and Ad, and between the strange and non-strange polarized sea densities 
are unnecessary and we have even considered the possibility that s(x) 7^ s(x) and 
As(x) 7^ As(x). 

Given the simplicity of the LO QCD analysis, we discuss where and when it is 
likely to be reliable, and stress the need to test the consistency of the LO treatment 
at each step. In this connection we have introduced the concept of a passive variable 
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in the experimentally measured observables. We have also suggested how one might 
estimate the errors induced in doing the LO analysis. 

In the NLO treatment we have shown how the expressions for the experimental 
observables can be much simplified by incorporating information from the reaction 
e^e~ -^ hX in a novel way. Determination of the polarized valence quark densities 
A-uv' and A(iy is shown to be relatively straight forward, as is the difference Am— A(i. 
However we argue that the determination of A-u, Arf, As + As and AG separately, 
from semi-inclusive DIS involving production of tt, K, A is unlikely to be successful, 
because of the limited range of Q^ available now and in the foreseeable future. It is 
suggested that the independent determination of AG from charm production is thus 
an essential element if A-u, Aii, As + As are to be determined accurately. Finally, 
motivated by the arguments that possibly s{x) ^ s{x) and As(a;) 7^ As(a;), we have 
demonstrated how, in principle, one can learn about s{x) — s{x) and As(a;) — As(a;). 

The procedure we have advocated poses a real challenge to the experimentalists, 
since it requires a control over the systematic errors involved in hadron detection 
efficiencies. The price paid in the current practice of considering certain ratios of 
cross sections in order to limit systematic errors, is enormous, and vast amounts of 
interesting and theoretically valuable information are thereby lost. We hope this 
paper will encourage efforts to proceed further. 
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